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Abstract 

Purpose Production of feed is an important contributor to 
life cycle greenhouse gas emissions, or carbon footprints 
(CFPs), of livestock products. Consequences of methodo¬ 
logical choices and data sensitivity on CFPs of feed ingre¬ 
dients were explored to improve comparison and 
interpretation of CFP studies. Methods and data for emis¬ 
sions from cultivation and processing, land use (LU), and 
land use change (LUC) were analyzed. 

Method For six ingredients (maize, wheat, palm kernel 
expeller, rapeseed meal, soybean meal, and beet pulp), 
CFPs resulting from a single change in methods and data 
were compared with a reference CFP, i.e., based on IPCC 
Tier 1 methods, and data from literature. 
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Results and discussion Results show that using more de¬ 
tailed methods to compute N 2 0 emissions from cultivation 
hardly affected reference CFPs, except for methods to de¬ 
termine NO^ leaching (contributing to indirect N 2 0 emis¬ 
sions) in which the influence is about -7 to +12 %. Overall, 
CFPs appeared most sensitive to changes in crop yield and 
applied synthetic fertilizer N. The inclusion of LULUC 
emissions can change CFPs considerably, i.e., up to 
877 %. The level of LUC emissions per feed ingredient 
highly depends on the method chosen, as well as on assump¬ 
tions on area of LUC, C stock levels (mainly aboveground C 
and soil C), and amortization period. 

Conclusions We concluded that variability in methods and 
data can significantly affect CFPs of feed ingredients and 
hence CFPs of livestock products. Transparency in methods 
and data is therefore required. For harmonization, focus 
should be on methods to calculate NO^ leaching and 
emissions from LULUC. It is important to consider LUC 
in CFP studies of food, feed, and bioenergy products. 

Keywords Carbon footprint • Feed ingredients - Feed 
production • Inventory data • Livestock products • Methods 

1 Introduction 

Environmental consequences of livestock production have 
received increasing attention over the last few years. Global 
warming, induced by emission of greenhouse gasses 
(GHGs), is one of the main problems addressed (Steinfeld 
et al. 2006). Livestock production contributes to global 
warming by emission of carbon dioxide (C0 2 ) from fossil 
fuel combustion and land use change (mainly deforestation), 
emission of methane (CH 4 ) from manure and enteric fer¬ 
mentation by ruminants, and emission of nitrous oxide 
(N 2 0) from manure storages and application of fertilizer 
for cultivation (Steinfeld et al. 2006; IPCC 2007; De Vries 
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and De Boer 2010). With livestock production being an 
important contributor to GHG emissions and the growing 
societal concern about global warming, GHG emissions 
from livestock production have become an imperative study 
object (Ellis et al. 2008; De Boer et al. 2011). 

Life cycle assessment (LCA) is an internationally accept¬ 
ed and standardized method (ISO 14043 2000) to evaluate 
GHG emissions of a product or production system. It eval¬ 
uates the use of natural resources and emission of pollutants 
along the entire life cycle of a product (Guinee et al. 2002; 
Rebitzer et al. 2004). Carbon footprint (CFP) assessment is a 
single-issue LCA focusing on emission of GHGs. 

The CFP of various livestock products has been calculat¬ 
ed, e.g., for milk (Haas et al. 2001; Thomassen et al. 2008, 
2009; Van der Werf et al. 2009; Flysjo et al. 2011a, b), pork 
(Basset-Mens and Van der Werf 2005), beef (Casey and 
Holden 2005; Beauchemin et al. 2011), chicken (Pelletier 
2008), and eggs (Mollenhorst et al. 2006; Dekker et al. 
2011). Such CFP assessments result in the identification of 
hotspots for GHG emissions along the production chain 
(Thomassen et al. 2009). A hotspot is a production stage 
with a high contribution to the environmental impact of a 
product. For most livestock products, this hotspot is feed 
production, including cultivation, processing, and transport 
stages. For milk, for example, production of feed explains 
around 45 % of the CFP (Thomassen et al. 2008; Van 
Middelaar et al. 2011); for pork, it is 60 %, and for chicken, 
even 80 % (Basset-Mens and Van der Werf 2005; Pelletier 
2008). Correct assessment of the CFP of feed ingredients, 
therefore, is an important aspect of CFP assessment of 
livestock products. 

To assess the CFP of feed ingredients, we need a harmo¬ 
nized method to calculate GHG emissions along the feed 
production chain. Variability in methods hampers compari¬ 
son of CFP results among studies (De Vries and De Boer 
2010). Particularly, accounting for emissions or C seques¬ 
tration from land use (LU) and land use change (LUC) 
appears complex. So far, there is no international consensus 
on a method to account for this, which increases variability 
in CFP studies (Prudencio da Silva et al. 2010; Cederberg et 
al. 2011; Flysjo et al. 2012). Exploring variability in meth¬ 
ods contributes to harmonization as it identifies the aspects 
that lead to differences between CFP studies. 

In addition to a harmonized method to calculate emis¬ 
sions, we need high-quality inventory data for each activity 
in the production chain, i.e., data on use of resources, 
emission of pollutants, and technical in- and outputs. Such 
inventory data can be subject to high uncertainty and vari¬ 
ability (Flysjo et al. 2011a). To improve LCA studies, in¬ 
sight into the relation between input data and the outcome of 
the study is required (Steen 1997; Sakai and Yokoyama 
2002). A sensitivity analysis shows for which data the 
outcome (e.g., the CFP of a product) is most sensitive. In 


other words, it shows which data should be considered first 
to improve the accuracy of an LCA study (Steen 1997). 

To improve comparison and interpretation of CFP studies, 
this study explored the effect of variability in methods and data 
sensitivity on CFPs of feed ingredients. We included emis¬ 
sions related to crop production and processing, and explored 
methods to account for GHG emissions or C sequestration 
from LULUC. Our objectives were: to give an overview of 
current methods that are used in CFP assessment of feed 
ingredients, to demonstrate consequences of methodological 
choices on final CFPs of feed ingredients, and to demonstrate 
sensitivity of CFPs of feed ingredients to technical in- and 
output data by performing a data sensitivity analysis. 

2 Material and methods 

2.1 Analysis framework 

Six feed ingredients were used to demonstrate consequences 
of methodological choices and data sensitivity on CFPs of 
feed ingredients, i.e., wheat, maize, soybean meal, palm 
kernel expeller, rapeseed meal, and beet pulp. We selected 
these ingredients because they are important ingredients in 
livestock concentrates, with major differences in nutritional 
value (Product Board Animal Feed 2008), and different 
production processes. 

To assess the CFP of feed ingredients, the following 
activities along the production chain are of importance: 
production of the system inputs (e.g., fertilizers, pesticides, 
energy resources), cultivation and harvesting of crop prod¬ 
ucts, drying and processing of crop products into single feed 
ingredients (this also includes the production of energy 
sources and auxiliary materials), processing of feed ingre¬ 
dients into a compound feed, and transport of unprocessed 
and processed products between all activities, up to the farm 
were the feed is used for livestock. This study included all 
activities up to the gate of the factory responsible for drying 
and processing of the single feed ingredients. 

The main GHGs emitted during production of feed ingre¬ 
dients are C0 2 , N 2 0, and, to a lesser degree, CH 4 (Duxbury 
1994). Production of system inputs, such as synthetic fertil¬ 
izers and energy resources, contributes mainly to C0 2 emis¬ 
sion, whereas N 2 0 emissions are most important in crop 
cultivation. Emission of CH 4 is minor and mainly related to 
peat soils (IPCC 2006). Emissions from LULUC are dom¬ 
inated by C0 2 . The CFPs of the six ingredients were com¬ 
puted by summing up emissions of these three gasses based 
on their equivalence factor in terms of C0 2 equivalents 
(100-year time horizon): 1 for C0 2 , 298 for N 2 0, and 25 
for CH 4 (IPCC 2007). 

Methods and inventory data for calculating CFPs of ingre¬ 
dients were collected from literature and by contacting 
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research institutes in France (INRA), Sweden (SIK), and the 
Netherlands (WUR). For wheat and rapeseed meal, data from 
several countries were used, resulting in nine data sets. 
Technical in- and output data for cultivation, drying, and 
processing of feed crops are included in the Electronic sup¬ 
plementary material 1, Table I. Yield and allocation factors of 
feed ingredients per feed crop are in the Electronic supple¬ 
mentary material 1, Table II. 

2.2 Reference CFP 

For each feed ingredient, a CFP in its most basic form was 
calculated, serving as a reference value to evaluate consequen¬ 
ces of methodological choices and data sensitivity. 
Computations of the reference CFPs were based on the fol¬ 
lowing assumptions. Emissions related to production of system 
inputs were based on life cycle inventories of the Ecoinvent 
database (2007). Production of seeds for sowing was not 
included. The amount of N from crop residues was based on 
IPCC (2006). Emissions of N 2 0 from crop cultivation were 
based on IPCC Tier 1 (IPCC 2006), which uses little or no 
country-specific data. Emissions related to drying and process¬ 
ing of ingredients were based on Ecoinvent (2007). For trans¬ 
port, an average distance per ingredient was used, based on 
country of origin (see Electronic supplementary material 1, 
Table I), whereas transport emission factors (EFs) were taken 
from Ecoinvent (2007). Emissions related to LULUC were not 
included in the reference CFPs. They were treated as a meth¬ 
odological choice. In case of a multiple output system, we used 
economic allocation. Economic allocation implies that the 
impact of a certain process is allocated to the various products 
based on their relative economic value. This type of allocation 
is mostly used in CFPs of feed products. Allocation factors are 
in the Electronic supplementary material 1, Table II. 

We demonstrated consequences of methodological 
choices and data sensitivity on CFPs of feed ingredients 
by comparing the CFP resulting from a single change with 
the reference CFP. Four categories of methodological 
choices were distinguished and are described in the follow¬ 
ing paragraphs: choices related to computation of emissions 
from cultivation (excl. LULUC), to emissions from LU, to 
emissions from LUC, and to emissions from processing. For 
emissions from cultivation and processing, we solely focused 
on emission calculations and not on the effect of changing the 
system boundaries, or allocation procedure. Although these 
aspects can have a large impact on the results, they have been 
subject to several other studies already (Flysjo et al. 2011b; 
Zehetmeier et al. 2012). For the data sensitivity analysis, the 
effect of a 10 % change in various inventory data on CFPs was 
examined, while keeping the other parameters constant. Data 
used for the data sensitivity analysis are in Table 1. The 
meaning and relevance of the data are described in the method 
sections below. 


2.3 Methods to compute GHG emissions from cultivation 
(excl. LULUC) 

To calculate GHG emissions from cultivation (excl. LULUC), 
we need methods to determine N 2 0 emissions from cultivation 
and to determine C0 2 emission from liming and urea fertiliza¬ 
tion. Emissions of N 2 0 from crop cultivation occur via a direct 
and an indirect pathway. Direct N 2 0 emission follows from 
microbial nitrification and denitrification of N in the soil. 
Indirect N 2 0 emissions involve N that is removed from soils 
via volatilization (e.g., ammonia (NH 3 ) or nitrogen oxide 
(NO*)), leaching, or runoff (e.g., nitrate (N0 3 ~)) (IPCC 2006). 
C0 2 from liming and urea fertilization occurs via dissolving of 
carbonates (C0 3 ) in C0 2 and water (H 2 0) (IPCC 2006). The 
literature review has revealed that such emissions are generally 
based on IPCC Tier 1, using an EF of 0.12 for limestone 
(CaC0 3 ), 0.13 for dolomite (CaMg(C0 3 ) 2 ), and 0.2 for urea, 
all expressed as kilograms of C0 2 -C per kilogram of product. 
No other methods for C0 2 emissions, therefore, were examined. 

2.3.1 N 2 0 emissions 

All peer-reviewed studies that calculated CFPs of feed 
ingredients that were found used IPCC (2006) to compute 
direct and indirect N 2 0 emissions. Some based their com¬ 
putations on general EFs as described in IPCC Tier 1, 
whereas others used country or fertilizer-specific EFs, or 
simulation models (Tier 2 and 3). These methods are pre¬ 
sented in Fig. 1. Direct N 2 0 emissions depend on the 
amount of inorganic nitrogen (N) available in the soil. In 
crop cultivation the available inorganic soil N increases due 
to the application of N fertilizers, the decomposition of crop 
residues, and the mineralization of soil N through LULUC. 
The latter is considered in the sections that describe the 
consequences of LULUC (§2.4 and §2.5). Only for the 
Netherlands, national inventory reports provide country- 
specific EFs to calculate direct N 2 0 emissions (in kilograms 
of N 2 0-N per kilogram N applied; Van der Hoek et al. 
2007), i.e., 0.005 for synthetic ammonium fertilizers, 0.02 
for manure (incorporating into the soil), and 0.01 for crop 
residue N (Van der Hoek et al. 2007). We therefore evalu¬ 
ated consequences of using these specific EFs on CFPs of 
Dutch feed ingredients only. Indirect N 2 0 emission is a 
function of volatilization of NH 3 and NO*, and leaching of 
NO^ (see Fig. 1). To compute volatilization of NH 3 and 
NO*, two other methods and for leaching ofNO^ three other 
methods were used besides Tier 1 (see Fig. 1). 

2.3.2 NH 3 and NO x volatilization 

We compared two methods to compute NH 3 and NO* volatil¬ 
ization with the reference situation. The reference situation 
uses one EF for synthetic fertilizer N and one for manure N 
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Table 1 Data used for data sensitivity analysis per category 

Cultivating and processing 

Land use 

Land use change 


Cultivating 

Processing 


Emissions per ha 

Emissions per feed ingredient 

Crop yield 

Synthetic fertilizer N 

Manure N 

Crop residues 

CaC0 3 (liming) 

Diesel use 

Emission factors 

Transport, feed crops 
Energy use drying 

Energy use processing 
Product yield 

Price (allocation factor) 

Reference soil C stock 

Soil C stock change factor 
Amortization period 

C stocks before LUC a 

C stocks after LUC a 

Soil N stock before LUC 

Soil N stock after LUC 

Amount of biomass burnt 

Total C stock change 

Area of LUC (ha) 
Amortization period 
Allocation factor b 


a Including C stocks in above- and belowground biomass, dead organic matter, and the soil 
b Allocation of LUC emissions to logging 


(IPCC Tier 1). The two other methods are: (1) using fixed 
fertilizer-specific EFs and (2) using a simulation model (see 
Fig. 1). For fertilizer-specific EFs for NH 3 volatilization per 
feed crop, see the Electronic supplementary material 1, 
Table III. The simulation model resulted in a country- 
specific, detailed prediction of NH 3 volatilization, taking into 
account the type of fertilizer, soil conditions, application tech¬ 
nique, and seasonal influences such as weather conditions 
(Karlsson and Rodhe 2002). This method was available for 
feed ingredients from Sweden only. In IPCC Tier 1, NO x 
volatilization is included in the EFs that are used; for the other 
two methods, NO x volatilization was based on Ecoinvent 
(2007), i.e., 0.21 multiplied by direct N 2 0-N emissions. 

2.3.3 NOT] leaching 

We compared three methods to compute NO^ leaching with 
the reference situation. The reference situation quantifies 
NO^ leaching as a fixed fraction of applied fertilizer N and 


crop residue N (IPCC Tier 1). The three other methods are: 
(1) the N field balance, (2) the NO^ leaching risk classes 
method, and (3) using a simulation model (see Fig. 1). The 
N field balance computes the difference between N inputs 
and N outputs at field level. This difference, also referred to 
as N surplus, is assumed to leach as NO^, although in 
practice, several other factors influence NO^ leaching. The 
NO^ leaching risk classes method assumes optimal fertil¬ 
ization and determines NCfT leaching by assigning crops to 
one of four leaching risk classes, based on type of crop, 
succeeding crop, duration of period without a crop, and 
postharvest soil N content (Basset-Mens et al. 2007). 
Postharvested soil N content is based on literature and 
expert’s opinion. Quantities of NCfT leaching per risk class 
are based on country-specific models. For France, for ex¬ 
ample, risk classes include 15, 40, 70, or 100 kg N0 3 ~-N/ha 
(Basset-Mens et al. 2007). This method was only available 
for feed ingredients from France and for soybean meal. The 
simulation model results in a country-specific, detailed 


Fig. 1 Methods for direct and 
indirect N 2 0 emissions, and for 
NH 3 +NO x volatilization and 
NO^ leaching (i.e., used to 
calculate indirect N 2 0 
emissions) in crop cultivation. 
a Used for NO^ leaching 
methods other than IPCC Tier 1 


N synthetic fert. 
N manure 
N crop residues 
N deposition a 
N fixation a 


NH 3 -N + NO x -N volatilization 

1. IPCC Tier 1 

0.1 * N syntfert + 0.2 * N manure 

2. Fertilizer specific EFs 

3. Simulation model 


Direct N 2 0-N emissions 

1. IPCC Tier 1 

0.01 * (TV syntfert + N manure + N cropres) 

2. Country-specific EFs 


Indirect N 2 0-N emissions 
1. IPCC Tier 1 

0.01 * (NH 3 + NO x ) + 0.0075 * N0 3 



N0 3 ' -N leaching 

1. IPCC Tier 1 

0.3 * (TV syntfert + N manure + N cropres) 

2. N field balance 

3. Leaching risk classes 

4. Simulation model 
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prediction of NO^ leaching, taking into account the type of 
fertilizer, soil conditions, application techniques, and 
ground water level (SEPA 2008). This method was only 
available for feed ingredients from Sweden. 

2.4 Methods to compute GHG emissions from land use 

In our study, LU refers to changes in management of crop¬ 
lands. LU can contribute to GHG emissions by affecting soil 
C stocks. An increase in soil C indicates removal of C0 2 
from the atmosphere (C sequestration), whereas a decrease 
indicates C0 2 emission. In addition, a decrease in soil C 
leads to N mineralization and hence N 2 0 emission (IPCC 
2006; Vellinga et al. 2004). It is assumed that when land use 
type and management system remain unchanged for deca¬ 
des, the soil C stock will no longer increase or decrease, and 
stabilization is reached. In our reference situation, therefore, 
no changes in soil C, and hence no LU emissions, were 
assumed. 

Parameters that affect soil C, and hence cause emissions 
from LU, are changes in the level of C inputs and changes in 
management practices that disturb the soil structure, such as 
tillage (Ogle et al. 2012). The level of C input highly depends 
on the amount of crop residues remaining on the field, 
depending on crop yield and crop residue removal. Manure 
application also is a source of C input. Changes in manage¬ 
ment practices that affect crop yield, such as a change in 
irrigation or fertilization regime, are related to LU emissions 
due to their indirect effect on the amount of crop residues 
(IPCC 2006). Changing to no tillage has been suggested as a 
strategy to decrease decomposition rates and increase C se¬ 
questration (Zotarelli et al. 2012). Recently, however, this 
effect has been questioned as a change in tillage system also 
can affect crop yield in a positive or negative direction, and 
hence C input as well (Ogle et al. 2012). 

Literature review shows that, so far, LU emissions have 
not been included in CLPs of feed ingredients. Assessing 
LU emissions requires detailed information on current and 
historical management practices, which is often not avail¬ 
able. furthermore, methods to calculate soil C stock 
changes have high levels of uncertainty. To gain insight into 
the potential consequences of including LU emissions on 
CLPs of feed ingredients, we calculated the effect of a 
change in tillage system on soil C stock levels based on 
IPCC Tier 1 methods and default values (IPCC 2006). To 
estimate the effect of a change in a tillage system, we need a 
reference tillage system. This was no tillage for soybeans 
and full tillage for all other feed crops. A change in tillage 
system, therefore, means changing to reduced or full tillage 
for soybeans and changing to reduced or no tillage for all 
other feed crops. Palm kernel expeller was excluded from 
the analysis as palm fruit is a perennial crop and does not 
require tillage. 


Emissions from LU depend on the level of soil C in the 
reference situation. Lor soybeans the reference soil C stock 
in the top layer (0-30 cm) was assumed to be 35 t C/ha; for 
all other crops, it was 32 t C/ha, based on estimates for C 
stocks in cropland after LUC from native vegetation with an 
average of 60 t C/ha (IPCC 2006). C stocks in soybeans 
were higher because as a default, no tillage is assumed to 
result in higher soil C stocks than full tillage (IPCC 2006). 
The C stocks were based on very rough estimates, but in line 
with 36 t C/ha in cropland in the Netherlands from Vellinga 
and Hoving (2011) and 28 t C/ha in cropland in Brazil from 
Cederberg et al. (2011). Emissions were amortized over a 
period of 20 years (IPCC 2006). 

2.5 Methods to compute GHG emissions from land use 
change 

In our study, LUC refers to transformation of noncropland, 
such as forest land, scrubland, and natural grassland, into 
cropland. LUC can contribute to GHG emissions by affect¬ 
ing C stocks in the ecosystem, including C stocks in above - 
and belowground biomass, dead organic matter, and soil 
organic matter. An increase in C stocks contributes to C 
sequestration, whereas a decrease contributes to C0 2 emis¬ 
sions. C0 2 emissions, for example, can occur from 
(incomplete) burning of aboveground biomass (e.g., defor¬ 
estation), from decay of biomass, and from changes in soil 
C. In addition, changes in soil C can lead to N 2 0 emissions 
(Vellinga et al. 2004), and burning of biomass leads to N 2 0 
and CH 4 emissions. 

To calculate LUC emissions related to crop cultivation, 
we need methods to estimate GHG emissions per hectare 
LUC (i.e., amount of emissions resulting from transforming 
1 ha of noncropland into cropland), and we need to allocate 
LUC to a specific crop (i.e., how many hectares are 
changed, which part of LUCs are allocated to which crop). 

2.5.1 Estimating GHG emissions per hectare of LUC 

Methods that estimate GHG emissions per hectare of LUC 
generally quantify changes in C stocks (Searchinger et al. 
2008; Leip et al. 2010; Cederberg et al. 2011). They vary in 
type of emissions accounted for and in time period over 
which changes in C stocks are examined (Electronic sup¬ 
plementary material 1, Table IV). We studied the conse¬ 
quences of including or omitting different types of 
emissions by evaluating the contribution of each type of 
emission to total LUC emissions per hectare, for situations 
relevant to feed crops. These were: changing tropical forest, 
scrubland, and natural grassland into annual cropland in 
Brazil (i.e., relevant for soybean meal) and changing tropi¬ 
cal forest into perennial cropland in Malaysia (i.e., relevant 
for palm kernel expeller). Calculations were based on IPCC 
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Tier 1 methods and default values. For C stocks in different 
land use categories, see Electronic supplementary material 
1, Table V. We assumed that part of the aboveground bio¬ 
mass was burned, i.e., 36 % of 160 t DM/ha biomass in 
tropical forest, 72 % of 14.3 t DM/ha biomass in scrubland, 
and 92 % of 5.2 t DM/ha biomass in grassland (IPCC 2006). 

We studied the consequences of a difference in time 
period over which changes in C stocks are examined by 
comparing the annual balance method (IPCC 2006), with 
the net committed emissions method (Fearnside 1997; 
Cederberg et al. 2011). The annual balance method is used 
most commonly and focuses on a specific time period, i.e., 
the moment that the land is cleared and used for another 
purpose, e.g., cropland. It does not include delayed emis¬ 
sions or C sequestration other than in the first year after 
LUC. The net committed emissions method encloses a 
longer time period and includes all delayed emissions and 
C sequestration that take place after the initial LUC. For 
soybean production in Brazil, for example, this method 
accounts for the fact that part of the land that was initially 
cleared for soybean (or pasture) production is abandoned 
after a few years. This abandoned land may regenerate into 
secondary forest, which can sequester C in biomass and soil, 
but it also means that more than 1 ha of land is changed to 
provide one new hectare of soybeans in permanent produc¬ 
tion. Differences between the two methods were analyzed 
for deforestation of tropical forest for soybean production in 
Brazil. Land use dynamics for the net committed emissions 
method were based on Macedo et al. (2012), assuming that 
15 % of the deforested land was abandoned in a later stage. 
As 34 % of the deforested land could not be classified into a 
land use category (Macedo et al. 2012), this 15 % was 
increased up to 23 %, assuming that similar transition prob¬ 
abilities hold for the unclassified category. We assumed that 
abandoned land regenerates into secondary forest and that 
the proportion of land that was initially deforested for crop 
production and later transformed into pasture was negligi¬ 
ble. Emissions per hectare of LUC were calculated similarly 
according to the annual balance method (IPCC 2006). 

2.5.2 Allocation of LUC to a specific crop 

To allocate LUC to a specific crop, we need to decide which 
crops are responsible for which part of the LUCs. Methods 
that are described in literature show high variation. We 
compared three methods. Method 1 focuses on direct LUC 
within a country or region and allocates emissions to the 
crops that are directly related to the LUC (Jungbluth et al. 
2007; Prudencio da Silva et al. 2010). Soybean meal and 
palm kernel expeller were the only two feed ingredients 
related to direct LUC. For soybean meal, we assumed that 
1 % of the soy produced in Central West Brazil comes from 
tropical forest, and 3.4 % comes from scrubland, whereas 


soy from South Brazil does not contribute to LUC 
(Prudencio da Silva et al. 2010). For palm kernel expeller, 
we assumed that 100 % of the palm area in Malaysia comes 
from tropical forest (Jungbluth et al. 2007). Calculation of 
emissions per hectare of LUC is described in the former 
section. For soybean meal, emissions per hectare were based 
on the annual balance method and the net committed emis¬ 
sions method, whereas for palm kernel expeller, the annual 
balance method was used only. Amortization period was 
20 years. 

As opposite to method 1, the following two methods also 
included indirect LUCs. Method 2 was based on Leip et al. 
(2010) and focuses on LUCs within a country, or country 
block (i.e., a group of countries), after which emissions were 
averaged for European Union (EU) countries and non-EU 
countries. For each country (block), the total area of LUC 
was determined for a specific time period, and emissions 
were allocated to the crops that showed an increase in total 
cropland in that time period, based on their relative contri¬ 
bution. Different types of LUCs were included, such as the 
transformation of natural grassland, scrubland, and tropical 
forest into cropland. LUC emissions were averaged, result¬ 
ing in one weighted value per crop (product) from the EU 
and one weighted value per crop (product) from non-EU 
countries. Emission calculations were based on the annual 
balance method but did cover the total area of LUC (i.e., all 
LUC was included). They were mainly based on IPCC 
default values (IPCC 2006) and for tropical forest and 
scrubland similar to method 1. Amortization period was 
20 years. 

Method 3 was based on Audsley et al. (2009). This method 
considers total LUC emissions worldwide, and allocates it to 
all agricultural land in use for commercial food production. 
Emissions were derived from Barker et al. (2007) and includ¬ 
ed GHG emissions and C sequestration from forestry only. C 
sequestration was included as the method accounts for affor¬ 
estation too. No amortization was applied. The method 
resulted in a single emission factor of 1.43 t C0 2 /ha of 
agricultural land. 

2.6 Processing of feed ingredients 

Feed ingredients can originate from crops directly (e.g., 
wheat and maize) or from industrial processing of crops 
(e.g., palm kernel expeller, soybean meal, rapeseed meal, 
and beet pulp). Ingredients that derive from industrial pro¬ 
cessing are often by-products from the biofuel or food 
industry. Rapeseed meal, for example, is a by-product from 
the processing of rapeseeds, whereas beet pulp is a by¬ 
product from the processing of sugar beets. To compute 
the CFP of a feed ingredient, therefore, methods that deal 
with multiple output systems are required. We used eco¬ 
nomic allocation in case of a multiple output system (see 
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Section 2.2). The impact of the processing stage, therefore, is 
determined by the amount and type of energy and auxiliary 
materials, and the emissions factors that are used. For alloca¬ 
tion, product yield after processing and price data are impor¬ 
tant. Overall, the main input for processing is energy 
(electricity, natural gas, diesel), facilitating processing stages 
as washing (sugar beets), crushing (oil seeds), and drying 
(grains, meal, and pulp). Hexane is often used in the oil 
industry as a solvent extraction, but generally the use of 
auxiliary materials is limited. Division of the processing stage 
into subprocesses increases the accuracy of CFP studies, but is 
often limited by lack of data (ISO 14043 2000). By changing 
the allocation factor of processing in the sensitivity analysis 
(see Table 1), we show for which feed ingredients subdivision 
is most important. A more detailed description of processing 
of feed ingredients can be found in Jungbluth et al. (2007). 

3 Results and discussion 

3.1 Reference carbon footprint 

Figure 2 shows the reference CFP of feed ingredients and 
the fractional contribution of different processes. For main 
products, such as maize and wheat, and for unprocessed 
products, N 2 0 emissions from cultivation and production 
of synthetic fertilizer N are by far most important contrib¬ 
utors (>65 %). For feed crops that use little synthetic fertil¬ 
izer N (e.g., legumes such as soybeans) or use different 
management practices that dominate emissions (e.g., high 
levels of irrigation such as palm fruit), this differs. Besides 
N 2 0 emissions and synthetic fertilizer N, production and 


combustion of diesel is quite important (10 %). Other 
aspects in cultivation (production of P 2 0 5 and K 2 0 fertil¬ 
izers, pesticides, and machinery) have a minor contribution 
only (<5 %). Emissions from drying and transport (about 
10 % of the CFPs) increase as the difference between DM 
content of harvested and dried product increases, or when 
transport distances increase. 

For by-products, such as palm kernel expeller and soybean 
meal, processing stages are important. Also for by-products, 
however, N 2 0 emissions from crop cultivation and production 
of synthetic fertilizer N are important contributors, except 
when the (economic) allocation factor for assigning emissions 
from crop cultivation to the feed ingredient is low and further 
processing is required (e.g., beet pulp). 

Only activities up to drying and processing of the single 
ingredients were included. For several ingredients, transport 
might have had a bigger impact, and CFP might have been 
higher, when all activities up to the country of final desti¬ 
nation would have been included. This counts especially for 
feed ingredients from tropical areas that are exported to 
Europe, such as soybean meal and palm kernel expeller. In 
Prudencio da Silva et al. (2010), for example, shipping of 
soybeans from South Brazil to the Netherlands contributed 
23 % to the total CFP of these soybeans. 

3.2 Cultivation and processing (excl. LULUC) 

3.2.1 Methods for N 2 0 emissions from crop cultivation 

Table 2 shows the reference CFP of feed ingredients and 
consequences of using more specific methods to calculate 
N 2 0 emissions from cultivation. Country-specific EFs for 
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Fig. 2 Reference CFP of feed ingredients and fractional contribution of different processes 
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Table 2 Reference CFP of feed ingredients (in kilograms C0 2 e/ton) and consequences of using more specific methods to calculate N 2 0 emissions 
from cultivation (in percent) 


Feed ingredient 

kg C0 2 e/t 

Direct N 2 0 emission 

NH 3 +NO x volatilization a 

N0 3 leaching a 





Country-spec. EFs (%) 

Fertilizer-spec. 
EFs (%) 

Simulation 
model (%) 

N field 
balance (%) 

Leaching risk 
classes (%) 

Simulation 
model (%) 

Maize (FR) 

507 


-2 


-5 

+6 


Wheat (FR) 

506 


-1 


-7 

-4 


Wheat (NL) 

502 

+2 

-2 


-3 



Wheat (SE) 

423 


-2 

-2 

-3 


-4 

Palm kernel expeller (MY) 

56 


0 


+9 



Rapeseed meal (FR) 

424 


-1 


-2 

-3 


Rapeseed meal (SE) 

405 


-1 

-2 

-2 


-3 

Soybean meal (BR) 

483 


0 


+ 1 

+12 


Beet pulp (NL) 

816 

0 

0 


0 




The reference CFPs are based on IPCC Tier 1 methods. No value means no data available 
EFs emission factors 
a Contributing to indirect N 2 0 emissions 


direct N 2 0 emissions were available for two feed ingre¬ 
dients only. Using this method changed CFPs with only 0 
and 2 %. For wheat (NL), a decrease in emissions resulting 
from a lower EF for synthetic fertilizer N compared to the 
reference was leveled out by an increase in emissions result¬ 
ing from a higher EF for manure (Van der Hoek et al. 2007). 
For feed ingredients that use another ratio of synthetic and 
organic fertilizer, the relative change could increase. 

Indirect N 2 0 emissions were computed as 0.1 x (NH 3 + N 
O x ) + 0.0075 x NCfy. This means that CFPs will be changed 
only when a change in method substantially alters the amount 
of NH 3 +NO x or NCfy. Using more specific methods to 
estimate volatilization of NH 3 +NO x changed CFPs only with 
2 %, whereas using more specific methods to quantify NOT 
leaching changed CFPs with -7 to +12 %. Based on these 
results, correct assessment of NCfy leaching is most important 
when calculating N 2 0 emissions from crop cultivation. 

3.2.2 Data sensitivity of inventory data 

Table 3 shows the reference CFP of feed ingredients and 
consequences of a 10 % change in inventory data and EFs. 
Overall, changing crop yield and synthetic fertilizer N 
changed CFPs most. Results correspond exactly with results 
in Fig. 2. This means that the higher the contribution of a 
certain aspect to the CFP of an ingredient, the higher the 
impact of the relative change. Crop yield is related to the 
contribution of cultivation in total. Therefore, CFPs of main 
products (maize and wheat), and (by-)products that have 
little emissions from processing (rapeseed meal), are more 
sensitive to a change in crop yield than products in which 
processing is more important (palm kernel expeller, soybean 


meal, and beet pulp). Similarly, we can explain the relative 
change of CFPs due to a 10 % change in amount of synthetic 
fertilizer N, affecting emissions from production and applica¬ 
tion, i.e., important contributors for most ingredients (see 
Fig. 2), and the relative change of the CFP of beet pulp due 
to a change in energy use for drying. Consequences of a 
change in product yield after processing are higher when the 
contribution of processes after processing are minor and when 
the product has a high allocation factor (e.g., soybean meal). 
Consequences of a 10 % change in price were highest for feed 
ingredients with a low allocation factor. Thus, for these ingre¬ 
dients, division of the processing stage into subprocesses can 
have a major impact on CFPs. A 10 % change in other 
technical in- and output data hardly affected CFPs. 
Regarding the EFs, results show that a 10 % change in the 
EF for direct N 2 0 emissions changed CFPs most (about 3 %). 

In cultivation, the quantitative order, thus relative impor¬ 
tance of inputs per crop type, does not vary between studies. 
For most feed ingredients, therefore, high-resolution data for 
crop yield and synthetic fertilizer N are most important for 
correct CFP assessment. For imported feed ingredients, 
means and distance of transport can be paramount. In gen¬ 
eral, higher accuracy in CFPs can be achieved by analyzing 
the relative contribution of different processes and validat¬ 
ing data for those processes that have a major contribution. 

3.3 Land use 

3.3.1 Methods for emissions from LU 

Table 4 shows the reference CFP of feed ingredients and 
consequences of a change in tillage system. To evaluate LU 
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Table 3 Reference CFP of feed ingredients (in kilograms C0 2 e/ton) and consequences of a 10 % change in inventory data and EFs (in percent) 


776 


Int J Life Cycle Assess (2013) 18:768-782 


bo 

a 

• iH 

xn 

in 

CD 

O 

o 

S-H 

a 

C 

t: 

o 

a 

CO 

§ 

H 


£ 

o 

• t-H 
> 

> 


=5 

U 


CD 

(N 

o 

u 

bD 


£ 

<D 

<D 

S-H 

W) 

£ 

• i-H 

CD 

CD 

Ph 


<D 

O 

S-H 

Ph 


o ^ 

2 « 
pp 


\ CO 

SB 8 

S IS 

(D ‘C 

C3 

W -3 O 


CO 

CP 

O 

Sh 

o 


t: 
o 

CP 
co 

2 CD 

& c£ 


I m SI) 

o .B 


pp 

w 


PC 

o 

c3 

CD 


+ cS 

£| 

2 * 
pp O N= 
£ £ € 


o 

CD 

Sh 


T3 

PP 

w 


o 

(D 

Sh 


Pp 

W 


(D 

co 

(D 


o 

CN 

£ 


o 

CN 

£ 


rp t>J) 

o .a 

o a 

ue 


CO 

(D 

2 

CP t 3 
O * 

lL CZ) 

U 2 


(D 

Sh —. 

3 ^ 
« 0 s - 
c3 

2 z 


.2 e 

I z 

Hu 

is) +; 


2 

73 

CP 

2 Co 
U O 


o o 


Cp pp CP CP O 

< < < < 

£ £ Z 


pp CP CP pp O 
£ £ £ £ 


pp 

£ 


o CP 
£ 


CN 


O O CP pp pp 

n I! I! 


oo O'n oo O'n t"- 


r- po 
o o 
m >r-) 


<N m 
O <N 
*/”) Xj- 


oi pi 

^ £h, 

N ^ 
C 3 jfn 

2 £ 


d S' 
£ SS 

+j +j 
c3 c3 
(D (D 
PP PP 

£ £ 


Sh 

2h 

73 

CP 

x 

(D 

73 


t-" r-" ■rf 


(N (N PO O 


(N O 


CP o O OP 
£ 


O O O CN OOOO 


o o o 


OOOOO OOOO 


o o 


m m m m <N mm 


(N O 


O CP <N CP 

< < 

£ £ 


<N O 


o o 


o o 


po po m po m po po 


oo oo t-- o 


po •'st m m po 

m <n o oo 

't 't d" X 


p4 

Pp 


c3 

(D 


“ « 
“ CQ 


c3 

(D 


03 

<D 


hJ 

£ 


(D 

P4 


o 

CD 

l CD 

CO 


IS| 

c3 a 
PP 


O 

(D 

(D 

co 

(D 

& 


Oh 

O 

CP 


SP 
c3 
(D 

O 

O O 

IS) PQ 


CD 

I 

o 

* 

u 

8- 
H—> 

o 

Ch 

Dh 

§ 


Springer 


methods, the default scenarios are used. Changing from full to 
reduced tillage changed CFPs by -1 to -15 %; changing from 
hill to no tillage by -1 to -28 % (see Table 4). Changes in 
CFPs were lowest for feed ingredients from France, because 
EFs varied with moisture regime and were relatively low for 
France, which has a dryer climate than the Netherlands and 
Sweden (IPCC 2006). Changing to reduced tillage resulted in 
C sequestration of about 90 kg C/ha/year for France and 
140 kg C/ha/year for the Netherlands and Sweden. When 
changing to no tillage, this was 220 kg C/ha/year for France 
and 260 kg C/ha/year for the Netherlands and Sweden. These 
numbers are in line with results found by Ogle et al. (2012). 
For soybean meal (using no tillage in the reference situation), 
changes were more pronounced, i.e., +55 % for changing to 
reduced tillage and +81 % for changing to full tillage, because 
EFs for tropical and wet climates were higher compared to 
European climates, and soybeans have a relative low yield 
compared to other feed crops. For soybeans, changing to 
reduced tillage resulted in a soil C loss of about 260 kg C/ha/ 
year, whereas for changing to full tillage, this was 385 kg C/ha J 
year. This is in line with results found by Zotarelli et al. (2012). 

Effects of a change in tillage system on soil C stocks have 
been questioned (Ogle et al. 2012). Ogle et al. (2012) showed 
that the final effect of a change in tillage system depends on a 
combination of crop type, climate, soil type, fertilization level, 
and other aspects and can vary between years due to variation 
in, e.g., weather conditions. In cold and wet climates, chang¬ 
ing from full to no tillage can even result in a decrease in soil C 
stock levels (Ogle et al. 2012). IPCC (2006) provides an 
uncertainty range along with their default values that displays 
this variation. Results shown here, therefore, are a first rough 
estimate of possible changes in CFPs when including FU 
emissions, but do not cover the complexity that is required 
for a detailed evaluation. For a detailed evaluation, all differ¬ 
ent aspects that influence crop yield, C input, and soil C stock 
levels need to be included and assessed. A change in crop 
yield is of particularly interest as this will also affect the 
allocation of emissions from cultivation. 

FU emissions due to a change in tillage system are 
nonrecurrent, whereas N 2 0 emissions from N application 
are annual. FU emissions were amortized over a period of 
20 years (IPCC 2006); thus, after 20 years, CFPs are no 
longer affected. This 20-year period is arbitrary and an 
estimation of the time that it takes to get to a new soil C 
balance. Including consequences of a change in tillage sys¬ 
tem means that this change is assumed to be permanent. If 
not, C0 2 that is sequestered from the atmosphere is emitted 
again as soon as the “old” tillage system is reimplemented. 

3.3.2 Data sensitivity in LU emissions 

To evaluate data sensitivity, results from a 10 % change in 
the default soil C stock change factor are compared. These 
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Table 4 Reference CFP of feed 
ingredients (in kilograms 

Ref. CFP (kg C0 2 e/t) 


Consequences of a change 

in tillage system (%) 


CCLe/ton) and consequences ot 
a change in tillage system (in 
percent) (default soil C stock 
change factor (-10 %, +10 %)) 

Feed ingredient 

Full tillage 

Reduced tillage 

No tillage 


Maize (FR) 

507 

-8 

(+7, -8) 

-19 

(-17, -20) 


Wheat (FR) 

506 

-8 

(+9, -9) 

-21 

(-19, -23) 


Wheat (NL) 

502 

-11 

(+4, -12) 

-20 

(-18, -23) 


Wheat (SE) 

423 

-15 

(+5, -16) 

-28 

(-25, -30) 


Rapeseed meal (FR) 

424 

-9 

(+10, -9) 

-21 

(-19, -23) 


Rapeseed meal (SE) 

405 

-12 

(+4, -13) 

-23 

(-20, -25) 


Beet pulp (NL) 

816 

-1 

(0, -l) 

-1 

(- 1 ,- 1 ) 



No tillage 

Reduced tillage 

Full tillage 



Soybean meal (BR) 

483 

+55 

(+49, +61) 

+81 

(+73, +89) 


results are presented between brackets (see Table 4). In the 
default scenario, the soil C stock is multiplied by the default 
stock change factor of 1.05 when changing from full to 
reduced tillage, in the case of feed crops from France 
(IPCC 2006). For feed crops from the Netherlands and 
Sweden, this factor is 1.08. A change of-10 % means that 
these stock change factors become <1, i.e., the soil C stock 
decreases, resulting in C0 2 emissions instead of C seques¬ 
tration. A -10 % change in the default stock change factor, 
therefore, increases CFPs, whereas a +10 % change 
decreases CFPs (see Table 4). As the stock change factor 
is subject to a lot of uncertainty (IPCC 2006), CFPs that 
include emissions from a change in tillage system should be 
interpreted carefully. For other feed ingredients, the default 
stock change factor was considerably higher than 1, and a 
10 % change did not change results from C uptake into C 
losses, or vice versa. In such cases, results from a 10 % 
change in stock change factor also apply to 10 % change in 
soil C stock level, whereas a 10 % change in amortization 
period resulted in slightly lower changes. This can be 
explained by the function (C stock change factor x C stock 
level/amortization period). Consequences increased with an 
increase in relative impact of LU emissions in CFPs (i.e., 
results at the default stock change factor) and therefore were 
highest for soybean meal. 

3.4 Land use change 

3.4.1 Methods for emissions from LUC 

Table 5 shows the contribution of different emissions to total 
LUC emissions per hectare for different land use transitions. 
Emissions from LUC are dominated by C0 2 emissions. 
When changing tropical forest or scrubland into cropland, 
the majority of the C0 2 emissions result from changes in 
aboveground C, belowground C, and soil C. When changing 
grassland into cropland, the majority of the C0 2 emissions 
result from changes in soil C, but C0 2 emissions from 


changes in above- and belowground C and N 2 0 from 
changes in soil N were still quite important. Excluding one 
of these emissions would result in underestimation of LUC 
emissions. When including emissions from burning, part of 
the C in biomass will be emitted as CO and CH 4 , and can 
therefore no longer be emitted as C0 2 . The net contribution 
of N 2 0 and CH 4 emissions from burning of biomass (i.e., 
after correction for foregone C0 2 emissions) was minor (see 
Table 5). Without this correction, emissions from deforesta¬ 
tion increased with about 10 t C0 2 e/ha, i.e., 1 % of total 
LUC emissions/hectare. 

There is little information on C stock levels in soils and 
belowground biomass in perennial croplands. The default 
soil C stock change factor for transformation of natural land 
into perennial cropland is 1, which means no change in the 
long term (IPCC Tier 1). This default value has a high 
uncertainty (50 %), and because in this case soil C losses 
are more likely than C sequestration, emissions from 
changes in soil C might be underestimated. In addition, 
C0 2 emissions from changes in belowground biomass C 
were not included, which is also expected to be an 
underestimation. 

Frequently, LUC emissions are amortized over a period 
of 20 years. When applying amortization, it seems correct to 
include land use transitions and C sequestration over the 
same period. For annual croplands, in which C sequestration 
is negligible, this will only affect results when part of the 
cropland is changed into another land use type or abandoned 
during the amortization period. This would mean that more 
than 1 ha of land is transformed to provide 1 ha of cropland. 
The net committed emissions method accounts for such land 
use transitions and delayed emissions and C sequestration 
after the LUC. When applying amortization, therefore, this 
method seems to be most suitable. Lack of information on 
land use transitions, however, can hamper its use. 

Using the net committed emissions method for changing 
tropical forest into annual cropland in Brazil resulted in an 
emission of 778 t C0 2 e/ha (775 t from a change in C stocks 
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Table 5 Total LUC emissions per hectare (in tons C0 2 e/heetare) and contribution of different emissions (in percent) for different land use 
transitions 

Country 

Original land use 

Final land use 

Total t 
C0 2 e/ha 

C0 2 from 
aboveground 
C (%) 

n 2 o+ch 4 

from burning 
biomass a (%) 

C0 2 from 
belowground 

c (%) 

C0 2 from 
dead organic 
matter C (%) 

C0 2 from 
soil C (%) 

N 2 0 from 
soil N (%) 

Brazil 

Tropical forest 

Annual cropland 

825 

63 

0 

23 

2 

11 

1 


Scmbland 

Annual cropland 

297 

47 

0 

19 

0 

31 

3 


Natural grassland 

Annual cropland 

128 

8 

0 

13 

0 

71 

7 

Malaysia 

Tropical forest 

Perennial cropland 

496 

97 

1 

NAY 

3 

0 

0 


NA V not available 

a Corrected for foregone C0 2 emissions from biomass due to emission of CO and CH 4 from burning 


and 3 t from burning including a correction for foregone C0 2 
emissions), compared to 825 t C0 2 e/ha for the annual balance 
method (see Table 5). Per hectare of permanent cropland, 
however, 1.23 ha is deforested. Hence, total LUC emission 
per hectare of permanent cropland is 957 t C0 2 e. 

Table 6 shows the reference CFP of feed ingredients and 
consequences of including LUC emissions using three 
methods. For method 1, the difference between the annual 
balance method and the net committed emissions method 
are given also. Method 1 focuses on direct LUC. Including 
direct LUC increased the CFP of soybean meal with 35- 
38 %, whereas the CFP of palm kernel expeller increased 
with 877 %. Method 2 and 3 also include indirect LUC. 
Including indirect LUC via method 2 (Leip et al. 2010) 
mainly affected the CFP of those ingredients that expanded 
their cultivation area over the last 10 years, i.e., rapeseed 
meal (change in CFP is about +70 %) and soybean meal 
(change in CFP is +632 %). Including indirect LUC in¬ 
creased especially the CFP of soybean meal, because meth¬ 
od 2 includes LUCs related to the expansion of soybean 
cultivation in 24 non-EU countries (blocks). This value, 
therefore, included the significant increases in area of 


soybean cultivation in the whole of Brazil, but also outside 
Brazil, for example in Venezuela (Leip et al. 2010). Method 
3 (Audsley et al. 2009) uses one single EF per hectare of 
land. The lower the yield per hectare and the higher the 
allocation factor, the higher the emissions per kilogram of 
feed ingredient. Including indirect LUC via method 3 
changed CFPs of wheat and maize with about +40 %. This 
is much more than with method 1 or 2, whereas for soybean 
meal and palm kernel expeller, the change in CFPs was less 
than for method 1 and 2. 

LUC emissions generally dominate CFPs, but the final 
change in CFPs varies between methods. There is no shared 
consensus, and the method chosen will greatly affect the 
outcome. The best method depends on the objective of the 
study. To encourage individual companies or countries to 
invest in sustainable production and to stimulate them to 
reduce deforestation, the method should focus on the direct 
link between products and LUC, i.e., method 1. Stimulating 
individual companies could lead to the combined demand of 
many actors for more sustainable production (i.e., no defor¬ 
estation) and hence to reduced deforestation in the long term 
(Weidema 2003). When the objective, however, is to 


Table 6 Reference CFP of feed 
ingredients (in kilograms 


CFP without FUC 

Consequences of including FUC emissions (%) 

C0 2 e/ton) and consequences ot 
including LUC emissions using 
three different methods, with for 
method 1 the difference between 
the annual balance method 
(AB) and the net committed 


kg C0 2 e/t 

Method 1 

Method 2 

Method 3 



AB NCE 



Maize (FR) 

507 

— — 

+3 

+31 

emissions method (NCE) (in 

Wheat (FR) 

506 

— — 

+8 

+40 

percent) 

Wheat (NF) 

502 

— — 

+8 

+31 


Wheat (SE) 

423 

— — 

+10 

+42 


Palm kernel meal (MY) 

56 

+877 NAV 

NAV 

+52 


rapeseed meal (FR) 

424 

— — 

+69 

+42 

Methods: 1= direct LUC (this 

rapeseed meal (SE) 

405 

— — 

+73 

+37 

study), 2=Leip et al. (2010), 3 = 

soybean meal (BR) 

483 

+35 +38 

+632 

+82 

Audsley et al. (2009) 

beet pulp (NF) 

816 

_ _ 

NAV 

+2 

NA V not available 
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emphasize that because of globalization of food and feed 
markets, the agricultural sector as a whole is responsible for 
deforestation, then the method should not differentiate be¬ 
tween direct and indirect LUC. In this case, every hectare of 
land used for commercial production purposes should be 
allocated a share of LUC emissions (method 3; Audsley et 
al. 2009). Method 3 will stimulate efficiency and increasing 
crop yield, and will favor feed crops from regions where the 
growth potential is highest due to optimal agro-ecological 
circumstances, because reducing land use requirements is 
the only option to reduce LUC emissions and hence CFPs. 
This method, however, does not provide a strong direct 
incentive to reduce deforestation. 

As LUC emissions dominate CFPs of feed ingredients, 
including these emissions might diminish the incentive to 
reduce emissions from cultivation and production other than 
from LUC. To avoid this, and because of high uncertainty 
and variation in calculating LUC emissions (the method is 
highly debated, and there is no shared consensus), emissions 
from LUC should be presented separately from other emis¬ 
sions in CFPs (Flysjo et al. 2012). Moreover, this seems 
correct because LUC emissions are nonrecurrent and only 
affect CFP for a certain period (i.e., dependent on amortiza¬ 
tion period), whereas other emissions from cultivation and 
processing recur annually. 

3.4.2 Data sensitivity in LUC emissions 

Table 7 shows the total LUC emissions per hectare for 
different land use transitions and consequences of a 10 % 
change in input data needed to calculate these emissions. 
Except for the amount of biomass burned, all inputs are 
directly related to the level of C stocks they refer to. This 
means that the higher the C stock, the more paramount the 
consequences of a 10 % change in this C stock will be. 
When changing tropical forest into cropland, for example, 
aboveground biomass in the initial land use was the largest 
C stock (see Electronic supplementary material 1, Table V), 


and a 10 % change in aboveground biomass, therefore, 
affected emissions most (see Table 7). Similarly, for scrub¬ 
land, this is the soil C stock in the initial land use, and for 
grassland, this is the soil C stock in both the initial and final 
land use. 

We examined the sensitivity of CFPs that include direct 
LUC emissions (i.e., method 1) to changes in input data 
necessary to calculate LUC emissions per feed ingredient. 
The more relevant the data for computing LUC emissions 
are (see Table 5), and the more dominant LUC emissions are 
in CFPs (see Table 6), the larger the effect of a change in 
data, with a maximum of 10 % (i.e., equal to the change in 
data). A 10 % change means that the aspect is highly 
relevant and that the CFP is completely determined by 
LUC emissions. As LUC emissions are almost completely 
determined by a change in C stocks (see Table 5), the 
consequence of a change in total C stock was comparable 
to the consequence of a change in the area of LUC (3 % for 
soybean meal and 9 % for palm kernel expeller in both 
cases). Consequences of a change in amortization period 
were slightly less, but about the same (3 % for soybean 
meal and 8 % for palm kernel expeller), which can be 
explained by the function (emissions per hectare x area of 
LUC/amortization period). A 10 % change in allocation 
factor means that 10 % of the aboveground biomass was 
allocated to logging, instead of no allocation. For perennial 
croplands, LUC emissions were for 97 % determined by a 
change in aboveground biomass C, whereas for annual 
cropland, this was 63 % (see Table 5). Therefore, a change 
in allocation factor will affect the CFP of a perennial crop 
more than the CFP of an annual crop. In this study, a change 
in allocation factor changed the CFP of soybean meal by 
1 %, whereas the CFP of palm kernel expeller was changed 
by 9 %. 

For accurate evaluation of LUC emissions, C0 2 from 
changes in aboveground biomass, belowground biomass, and 
soil C should be included. For correct interpretations and 
comparisons of LUC emissions per feed ingredient, it is 


Table 7 Total LUC emissions for different land use transitions (in tons C0 2 e/heetare) and consequences of a 10 % change in different input data (in 
percent) 


Country Initial land use Final land use Total Aboveground Biomass Belowground Dead organic Soil C stocks Soil N stocks 




emissions 

biomass 

C stocks 

burnt 

biomass 

C stocks 

matter C stocks 







t C0 2 e/ha 

Initial 

(%) 

Final 

(%) 

Initial 

(%) 

Initial (%) 

Initial (%) 

Initial 

(%) 

Final 

(%) 

Initial 

(%) 

Final 

(%) 

Brazil Tropical forest 

Annual cropland 

825 

6 

NAP 

0 

2 

0 

3 

1 

0 

0 

Shrub land 

Annual cropland 

297 

5 

NAP 

0 

2 

0 

8 

5 

1 

0 

Natural grassland Annual cropland 

128 

1 

NAP 

0 

1 

NAP 

17 

8 

2 

1 

Malaysia Tropical forest 

Perennial cropland 

496 

10 

1 

0 

NAV 

0 

4 

4 

1 

1 


NAP not applicable, NAV not available 
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equally important to consider assumptions on the area of LUC, 
as C stock levels and amortization period. Assumptions about 
logging can be important too, especially for perennial 
croplands. 

3.5 General discussion 

We did a sensitivity analysis to identify for which data the 
outcome of the CFP studies are most sensitive. Such infor¬ 
mation can improve the accuracy of CFP studies as it points 
out which data should be considered first. The effect of a 
change in input data on the outcome of the study is deter¬ 
mined by the magnitude of the change. Changing the input 
data with 50 %, or 5 %, instead of 10 %, however, does not 
change the priority of the input data and will therefore not 
affect the conclusions of the sensitivity analysis. 

Results do not give insight into the effect of data uncer¬ 
tainty, which refers to uncertainty due to inaccurate meas¬ 
urements, or lack of data, or data variability, which refers to 
variation in the real world, e.g., temporal and spatial varia¬ 
tion (Huijbregts 1998). An uncertainty analysis requires 
information on distribution and data quality indicators, and 
can be performed with, for example, a Monte Carlo analysis 
(Heijungs and Huijbregts 2004). The effect of data variabil¬ 
ity can be large, but is not a matter of lack of data quality or 
knowledge. Crop yield, for example, varies greatly between 
countries, but also within countries and between years 
(FAOSTAT 2010), and has a large impact on the CFP of 
feed ingredients. The same accounts for application of ma¬ 
nure and fertilizers, including limestone. The effect of data 
uncertainty and variability, therefore, can be much larger 
than results shown by our sensitivity analysis. It is important 
to realize this when comparing CFP studies. Particularly, 
emissions from LULUC can vary greatly due to high levels 
of data uncertainty. 

4 Conclusions 

We explored the consequences of methodological choices 
and data sensitivity on CFPs of feed ingredients for emis¬ 
sions from cultivation and processing, LU, and LUC. 
Calculation methods for direct and indirect N 2 0 emissions 
from cultivation were consistent among studies, whereas 
differences in methods to calculate NH 3 and NO x volatili¬ 
zation (contributing to indirect N 2 0 emissions) hardly af¬ 
fected CFPs. Differences in methods to calculate N0 3 _ 
leaching (also contributing to indirect N 2 0 emissions), how¬ 
ever, can affect CFPs considerably. High-resolution data 
were most important for crop yields and the quantity as well 
as the type of synthetic nitrogen fertilizer. For by-products, 
data on processing and transport can be paramount. Higher 
accuracy in CFPs can be achieved by analyzing the relative 


contribution of different processes and validating data for 
the most important parameters, e.g., yield and N fertilizer 
data. 

We explored the consequences of including LU emis¬ 
sions (i.e., emissions due to a change in management prac¬ 
tices) by assessing the effect of a change in tillage system. 
Results show that changing to no tillage can potentially 
reduce CFPs. For a detailed evaluation, however, all aspects 
that affect crop yield, C input, and soil C stock levels should 
be included. For accurate evaluation of LUC emissions, 
C0 2 from changes in aboveground biomass, belowground 
biomass, and soil C should be included. The net committed 
emissions method seems to be most appropriate when ap¬ 
plying amortization: C stock changes and land use transi¬ 
tions are accounted for preferably over the same period as 
the amortization period. For allocating LUC to different 
crops, the objective of the study is important, and the meth¬ 
od will greatly affect results. LULUC emissions should be 
presented separately from other emissions, because there is 
no consensus about the method to calculate these emissions, 
and LULUC emissions are nonrecurrent, whereas other 
emissions reoccur annually. To compare LUC emissions 
per feed ingredient, the area of LUC, C stock levels, and 
amortization period should be considered. Assumptions 
about logging can be important too, especially for perennial 
croplands. 

Variability in methods and data can considerably affect 
CFPs of feed ingredients and hence CFPs of livestock 
products. Transparency in methods and data are necessary 
to distinguish between actual differences and differences 
caused by methods and data. For harmonization, focus 
should be on methods to calculate N0 3 ~ leaching and emis¬ 
sions from LULUC. It is important to consider LUC in CFP 
studies of food, feed, and bioenergy products. 
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